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ABSTRACT 

Context. We contribute to developing dynamical models of the motions of Uranus’ main satellites. 

Aims. We determine accurate positions of the main satellites of Uranus: Miranda, Ariel, Umbriel, Titania, and Oheron. Positions of 
Uranus, as derived from those of these satellites, are also determined. The observational period spans from 1992 to 2011. All runs 
were made at the Pico dos Dias Observatory, Brazil. 

Methods. We used the software called Platform for Reduction of Astronomical Images Automatically (PRAIA) to perform a digital 
coronography to minimise the influence of the scattered light of Uranus on the astrometric measurements and to determine accurate 
positions of the main satellites. The positions of Uranus were then indirectly determined by computing the mean differences between 
the observed and ephemeris positions of these satellites. A series of numerical filters was applied to filter out spurious data. These 
filters are mostly based on (a) the comparison between the positions of Oheron with those of the other satellites and on (b) the offsets 
as given by the differences between the observed and ephemeris positions of all satellites. 

Results. We have, for the overall offsets of the five satellites, -29 mas (±63 mas) in right ascension and -27 mas (±46 mas) in 
declination. For the overall difference between the offsets of Oberon and those of the other satellites, we have ±3 mas (±30 mas) in 
right ascension and -2 mas (±28 mas) in declination. Ephemeris positions for the satellites were determined from DE432±ural 11. 
Comparisons using other modem ephemerides for the solar system - INPOP13c - and for the motion of the satellites - NOE-7-2013 
- were also made. They confirm that the largest contribution to the offsets we find comes from the motion of the barycenter of the 
Uranus system around the barycenter of the solar system, as given by the planetary ephemerides. For the period from 1992 to 2011, 
our final catalogues contain 584 observed positions of Miranda, 1,710 of Ariel, 1,987 of Umbriel, 2,588 of Titania, 2,928 of Oberon, 
and 3,516 of Uranus. 

Key words. Astrometry - Methods: data analysis - Catalogs - Planets and satellites: individual: Uranus, Miranda, Ariel, Umbriel, 
Titania, Oberon 


1. Introduction 

Accurate positions from gound-based CCD images of small bod¬ 
ies and natural satellites in the solar system are an important 
tool for developing accurate dynamical models of their motions 
(Lainey 2008) and to support (see Desmars et al. 2015, for in¬ 
stance) further investigations of these bodies and their surround¬ 
ing environments by other observational techniques and meth¬ 
ods, such as stellar occultations (Braga-Ribas et al. 2014) In 


* Based on observations made at the Pico dos Dias Observatory/LNA. 
** Tables with the positions of Uranus and its satellites, as well 
as with the X and Y CCD coordinates of the observed satellites 
and reference stars, are available at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz- 
bin/qcat?J/A±A/???/A??? 


addition, the knowledge of the dynamics of natural satellites is 
important to constrain models of formation and evolution of the 
solar system (Nogueira et al. 2011) , to understand their internal 
structure (Lainey et al. 2009), and to more general studies of the 
physics of a particular planetary system (see Jacobson 2014, for 
the Uranus system) 

A summary of ground-based techniques used in the astrome¬ 
try of solar system bodies, along with their respective accuracies, 
can be found in Arlot (2008) Stellar occultations can also be in¬ 
cluded among those techniques since they are also a source of 
accurate positions (see Person et al. 2006; Sicardy et al. 2006; 
Widemann et al. 2009; Assafin et al. 2010; Benedetti-Rossi et al. 
2014) In the case of Pluto, for instance, one could take advan¬ 
tage of the fact that its photocenter displacement due to the pres- 
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Table 1. Characteristics of the telescopes 


Table 2. CCD detectors 


Telescope 

f 

(mm) 

f/ 

scale 

("/mm) 

System 

Perkin-Elmer 

15752 

f/10 

13.09 

RC 

Boller&Chivens 

8222 

f/13.49 

25.09 

RC 

Zeiss 

7500 

f/12.5 

27.50 

Cass 


Columns are manufacturer, focal length, f-number, pixel scale, 
and telescope system (Ritchey-Chretien - RC or Cassegrain - 
Cass). 


ence of Charon is not of concern for the astrometry from stellar 
occultations. 

Here, we aim at contributing to dynamical models of the mo¬ 
tions of the main satellites of Uranus: Miranda, Ariel, Umbriel, 
Titania, and Oberon. This contribution is made by determining 
accurate positions determined from CCD observations carried 
out at the Pico dos Dias Observatory in the period 1992-2011. 

Uranus was not directly observed in this work. However, po¬ 
sitions of Uranus, indirectly determined from those of the satel¬ 
lites, are also provided since they are relevant to improve its or¬ 
bit. In fact, orbits in modern planetary ephemerides are not uni¬ 
formly known. Those of Uranus, Neptune, and Pluto have uncer¬ 
tainties of several thousand kilometers, two orders of magnitude 
greater than those for Jupiter and Saturn. The situation is even 
more dramatic when the comparison is made with the inner so¬ 
lar system, where orbits are known to subkilometric accuracy 
(see Folkner et al. 2014, for further details) 

This work represents an extension of that by Veiga et al. 
(2003) (hereafter V03) and follows the efforts of many others 
that determined optical positions within the system of Uranus 
from the ground. In addition to traditional astrometry, a no¬ 
ticeable contribution comes also from the observation of mu¬ 
tual events, where Hidas et al. (2008) were the first to report 
such an analysis to our knowledge that involved the satellites 
of Uranus. Comprehensive lists of observational sources for the 
system of Uranus can be found, for instance, in Jacobson (2014) 
and |Emelyanov & Nikonchuk] ( |2013] l. 

Next, in Sect. 2, we describe the site and instruments. 
Data reduction procedure and the results and analysis are pre¬ 
sented in Sects. 3 and 4, respectively. Comparisons with other 
ephemerides are given in Sect. 5. Section 6 argues the extension 
of the astrometry of this work to older epochs using results from 
V03, and Sect. 7 presents our position catalogues. Comments 
and conclusions are given in Sect. 8. 


Type 

Size (pixels) 

Size (microns) 

A (nm) 

(a) 

2048 X 4608 

13.5 X 13.5 

500 

(b) 

2048 X 2048 

13.5 X 13.5 

600 

(c) 

1024 X 1024 

24 x24 

650 

(d) 

2048 X 2048 

13.5 X 13.5 

491 

(e) 

1024 X 1024 

13.0 X 13.0 

560 


CCD types are (a) Marconi CCD42-90-0-941, (b) Marconi 
CCD42-40-1-368, (c) SITe SI003AB, (d) Andor iKon-L, and (e) 
Andor iXon™. The other columns give the CCD dimensions, 
pixel sizes, and the wavelength of the CCD maximum quantum 
efficiency. 
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> > 
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« 
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Fig. 1. Image containing Uranus and its main satellites, as obtained on 
18 July, 1992, with the Perkin-Elmer telescope. Upper left (right) panel: 
inverted colour map image before (after) the application of the digital 
coronography procedure. Lower left (right) panel: respective Miranda 
flux profile along a segment placed on the line joining the centres of 
Uranus and Miranda. This segment is the same in the two lower panels. 
Image orientation is north up and east left. Exposure time is 9s. All other 
objects in the image that do not belong to the system of Uranus are field 
stars. 


2. Site and instruments 

All observations were made with three telescopes located at the 
Pico dos Dias Observatory (lAU code 874, A = -45°34'57" W, 
(j) - -22'’32'04", and h=1810.7m), run by the Laboratorio 
Nacional de Astroffsica/MCTQ Their characteristics are sum¬ 
marised in Table[2 All of them have equatorial mounts. 

A number of detectors, as given by Table|^ were installed at 
the telescopes presented in Table [T] Typically, our images from 
the Perkin-Elmer telescope have a field of view (EOV) of 3' to 
6', and for the other telescopes, the EOV is 5' to 11'. 

Eilters B, V, R, and /, compatible with the Johnsons-Cousins 
system, were used. Most of our observations, however, were 
made without filter. 

* http://www.lna.br/. Page in Portuguese. 
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3. Data reduction 

All raw images were first corrected for bias and flatfield with 
IRAE (Tody 1993) Then, a digital coronography procedure (see 
Assafin et al. 2008, 2009) was applied to minimise the effect of 
the scattered light from Uranus on the positions of the satellites. 

Eigure[T]illustrates the performance of this correction. In this 
figure, the angular separation between the centres of Uranus and 
Miranda is 9.5". Although the centroid determination of any ob¬ 
ject contaminated by the scattered light of Uranus profits from 
the digital coronography, this procedure is evidently more rele¬ 
vant to Miranda in our images. 

The final step, that is, the determination of the positions of 
the satellites of Uranus in all images, was fully implemented 
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Fig. 2. Distribution of the number of observations per year. 


with the Platform for Reduction of Astronomical Images Au¬ 
tomatically (PRAIA) packag^ (Assafin et al. 2011) 

The results presented here were obtained from the analysis 
of a total of 4,532 images, all telescopes and detectors shown 
in Tables and 1^ included, distributed along 138 nights in the 
period 1992-2011 (see also Fig.|^. These counts exclude images 
(23% of the total) where positions, for whatever reason, were 
not obtained. We used the Fourth US Naval Observatory CCD 
Astrograph Catalog (UCAC4 - Zacharias et al. 2013) as the 
practical representative of the International Celestial Reference 
System (ICRS - Arias et al. 1995) 

Figure]^ also shows the absence of observations in 2000 and 
2007. In 2000, we had no observational runs at the Pico dos Dias 
Observatory. In 2007, all observations (more than 18,000) of the 
Uranus system were aimed at determining dynamical parame¬ 
ters from light curves obtained during mutual events among the 
satellites (see Assafin et al. 2009) As a result of their short expo¬ 
sure times (1-3 seconds) and small FOV (usually 2'x2' or less), 
they were not used in this work. 


4. Results and data analysis 

The observed positions of the satellites were first compared 
to those provided by ephemerides from the JPL; the planetary 
ephemeris DE432 and ural 1 1[^ 

By means of the SPICE system toolkits (Acton 1996), 
DE432 was used to link the geocenter to the barycenter of 
the Uranus system, whereas ural 11 was used to link this lat¬ 
ter to the centre of Uranus and those of its main satellites. The 
link between the geocenter and the topocenter was obtained 
with subroutines from the Standards of Fundamental Astron- 
om}lj(SOF.M and from the Naval Observatory Vector Astrome¬ 
try Softwar^(NOVAS). The sum of these links (or vectors) pro¬ 
duced the astrometric ephemeris positions of Uranus and of its 
main satellites for an observer at the Pico dos Dias Observatory. 

The differences in the sense observed minus ephemeris posi¬ 
tions are referred to as offsets in position, following the nomen¬ 
clature commonly used in astrometric works. These offsets were 
determined and used to filter out spurious data. Details on this 

^ The digital coronography procedure is also part of this package. 

^ both at http://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/ 
http://www.iausofa.org 

^ http://aa.usno.navy.mil/software/novas/novas_info.php 


filtering are given below. We note that, except for item 1, no fur¬ 
ther consideration was made with respect to the telescope, filter, 
or detector used. 


1. All satellite positions whose centres were closer than 8" (for 
observations made with the Perkin-Elmer telescope) or 7" 
(for observations made with the Boller&Chivens or Zeiss 
telescopes) to the centre of Uranus were discarded. This was 
done to avoid cases where the digital coronography might 
fail to produce reliable results. 

2. All images with astrometry derived from less than five refer¬ 
ence (UCAC4) stars were discarded. This was done avoid a 
poor representation of the celestial frame. 

3. Eor each night and each satellite, an iterative cr-clipping pro¬ 
cedure filtered out all those offsets that differed from the 
mean offset in right ascension or declination by more than 

2- cr; this way, we eliminated poor measurements. 

4. To all observations of a given satellite, an iterative cr-clipping 
procedure filtered out all those residuals in right ascension 
and declination, as given by the difference between the off¬ 
sets of the given satellite and those of Oberon in the same re¬ 
spective images that differed from the mean residual by more 
than 3-cr. This was a coherence check, assuming that the or¬ 
bits of the satellites around the barycenter of the Uranus sys¬ 
tem are well known. This is discussed below. Oberon was 
chosen as reference because its positions were determined in 
all the 4,532 images mentioned earlier and, most relevant, it 
is the most distant and the slowest among the satellites stud¬ 
ied here (as a consequence, Oberon has the most realiable 
orbit of the satellites studied here). 

5. To all observations of a given satellite, an iterative cr-clipping 
procedure filtered out all those offsets in right ascension and 
declination that differed from the mean offset by more than 

3- cr. We did this to use the information from the whole data 
set to filter out poor representations of the celestial frame that 
survived item 2 and that can only be identified when a large 
portion of the sky is considered. 


These filters were applied one after the other in the sequence 
they have just been presented. The most severe of them, that is, 
the one that eliminated the most, was that described in item 3. 
The least severe procedure was that described in item 5. 

Eigure|^shows the number of observations per year after ap¬ 
plying all filters. Erom comparing Eigs.[^and[^ it is clear that 
positions derived from observations made after 2000 were more 
affected by the application of the filters than those made ear¬ 
lier. The reason for this can be easily explained by Eig. After 
2000, the planes of the orbits of the satellites were approaching 
an edge-on view so that observations became more complicated, 
and in addition, more satellite positions were filtered out because 
of their proximity to Uranus. 

Tables]^ to 1^ show yearly offsets for each satellite. Table ^ 
clearly shows the decrease in the number of positions of Miranda 
as the plane of its orbit approaches an edge-on view. Eigures ^ 
to 1^ are provided to illustrate these tables by showing all final 
offsets for each satellite as a function of time. We use an asterisk 
to indicate multiplication by the factor coscJ throughout the text. 

The standard deviations, as shown in Cols. 4 and 5 in Ta¬ 
bles to peak in right ascension. Although it is not easy to 
clearly identify the reason for this feature, some hypothesis may 
be investigated. 


® Generated by the Uranus Viewer Tool, PDS Rings Node, http://pds- 
rings.seti.org/ 
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1995 2000 2005 2010 

Date (year) 

Fig. 3. Distribution of the number of observations per year after apply¬ 
ing filtres 1 to 5. 

Table 3. Results for Miranda 


Year 

Aa* Ai5 
(mas) 

o~ s 

(mas) 

ea* es 
(mas) 

# 

positions 

1992 

33 

1 

66 

50 

53 

46 

54 

1993 

-32 

-22 

88 

60 

48 

54 

71 

1994 

60 

-41 

57 

52 

50 

51 

54 

1995 

16 

-2 

56 

62 

45 

48 

158 

1996 

-74 

-3 

93 

58 

59 

60 

116 

1997 

-89 

-22 

87 

40 

62 

44 

36 

1998 

14 

5 

38 

49 

49 

33 

24 

1999 

-116 

35 

62 

59 

44 

51 

22 

2001 

38 

-69 

48 

31 

30 

25 

19 

2002 

-63 

23 

149 

58 

53 

53 

6 

2004 

-271 

-83 

— 

- 

52 

34 

1 

2006 

-289 

17 

5 

23 

57 

60 

3 

2009 

-103 

25 



60 

66 

1 

2011 

-164 

38 

107 

61 

47 

49 

19 


Columns: year of observation; mean offset for the given year in 
right ascension; mean offset for the given year in declination; 
standard deviation of the offset in right ascension; standard devi¬ 
ation of the offset in declination; mean value of the standard de¬ 
viations in right ascension for the given year, as derived from the 
differences observed minus catalogue positions for the reference 
stars; mean value of the standard deviations in declination for 
the given year, as derived from the differences observed minus 
catalogue positions for the reference stars; number of positions 
that yielded the results shown in Cols. 1 to 7. All angular units 
are in mas (milliarcsecond). 





TimelUTC): 1999-JUL-01 



Time (UTC): 2003-JUL-01 



Time(UTC): 2007-JUL-01 



Time (UTC): 2011-JUL-01 


Fig. 4. Evolution of the pole position of Uranus, with the vie\4^ of the 
satellite’s orbit planes gradually changing from close to face-on (1992) 
to edge-on (2007). 


Known small mechanical problems that affect the telescope 
tracking system might be responsible for a degradation in the 
quality of the measurements in right ascension. However, in this 
case, standard deviations in right ascension should be systemat¬ 
ically larger than those in declination, and this is not verified. 

Another source of larger uncertainties are occurrences of in¬ 
correct timing inserted in the image headers. This mostly affects 
the right ascension measurements and is known to have hap¬ 
pened in the past (before 2000). Figure [T0| illustrates this point. 
We note that both mechanical and timing problems may slightly 
increase or decrease the values in right ascension in our case. We 
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therefore expect on average larger standard deviations and not a 
systematic effect. 

The only points above the threshold of 55 mas in Fig. [T^ 
are those in right ascension and not far from opposition, where 
the movement of Uranus (most of it along right ascension) is 
faster. The value of 55 mas is arbitrary to some extent and was 
chosen because the standard deviations in declination (Col. 5 in 
Tables [^to0 are almost all smaller than it. Therefore, some tim¬ 
ing problem is a plausible and probable explanation. Figure 
also shows a number of observations far from opposition that 
are therefore helpful for determining the heliocentric distance of 
Uranus and its satellites. 
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Miranda JPL 



Date (year) 
Miranda JPL 



Ariel JPL 



Ariel JPL 



Date (year) 


Fig. 5. Distribution of the offsets for Miranda as a function of time. 


Fig. 6. Distribution of the offsets for Ariel as a function of time. 


Table 4. Results for Ariel 


Year 

Aa* A6 
(mas) 

(mas) 

ea* es 
(mas) 

# 

positions 

1992 

8 

-13 

34 

20 

54 

48 

75 

1993 

-12 

-28 

19 

20 

49 

53 

95 

1994 

36 

-35 

40 

43 

50 

47 

56 

1995 

28 

-23 

43 

38 

50 

50 

206 

1996 

-40 

-4 

74 

43 

54 

51 

245 

1997 

-37 

-25 

89 

27 

56 

45 

78 

1998 

-4 

4 

37 

17 

50 

36 

28 

1999 

-69 

-8 

77 

53 

48 

51 

33 

2001 

-14 

-30 

39 

38 

44 

43 

46 

2002 

2 

-6 

32 

30 

52 

49 

91 

2003 

-16 

-25 

82 

69 

65 

68 

35 

2004 

-22 

-31 

65 

51 

57 

55 

102 

2006 

-79 

40 

32 

8 

61 

46 

2 

2010 

-168 

40 

27 

7 

79 

58 

15 

2011 

-61 

-30 

48 

38 

62 

58 

603 


Columns: same as Tabled 


Another interesting point about the offsets is shown in 
Fig. where it is possible to note that, except for Miranda, the 
yearly offsets in right ascension and declination are similiar. The 
overall offsets for all five satellites are shown in Table [3 This is 
an indication that the offsets we find are mostly due to the plan¬ 
etary ephemeris DE432, that is, due to the ephemeris position of 
the barycenter of the Uranus system. Support for this result is 
also obtained from Cols. 2 and 3 in Table |9] 

At this point, it is interesting to check whether our observed 
positions might also be an improvement to models of the satellite 
motions. A simple way of checking this is given by Figs.[T2|to 


16 In them, each symbol is the mean of 100 offsets, and the 


length of the error bars, above and below the respective dots or 
triangles, is three times the error of the mean. 


Table 5. Results for Umbriel 


Year 

Aa* Ad 
(mas) 

fa* fij 

(mas) 

Ca* eg 
(mas) 

# 

positions 

1992 

8 

-13 

31 

20 

54 

47 

60 

1993 

-11 

-30 

21 

25 

49 

54 

115 

1994 

36 

-54 

37 

48 

51 

49 

68 

1995 

27 

-24 

39 

45 

49 

49 

230 

1996 

-36 

-7 

66 

48 

56 

54 

259 

1997 

-22 

-23 

77 

20 

56 

47 

68 

1998 

-8 

15 

32 

15 

52 

36 

23 

1999 

-63 

-3 

73 

57 

49 

51 

36 

2001 

-18 

-33 

36 

35 

36 

37 

60 

2002 

-9 

-6 

36 

37 

54 

51 

122 

2003 

-29 

-13 

72 

43 

62 

58 

91 

2004 

-31 

-29 

55 

43 

59 

53 

no 

2005 

-18 

-69 

11 

10 

43 

54 

19 

2006 

-45 

-38 

40 

49 

56 

58 

61 

2009 

17 

-118 

22 

30 

61 

61 

30 

2010 

-117 

30 

36 

12 

72 

66 

99 

2011 

-54 

-47 

56 

48 

57 

52 

536 


Columns: same as Table|3] 


It is possible to note systematic deviations in the offsets of 
both right ascension and declination as a function of the true 
anomaly. Positions from observations made with the satellite or¬ 
bits close to edge-on tend to be more affected by systematic ef¬ 
fects than those from observations made with the satellite orbits 
close to face-on. However, both sets of positions present system¬ 
atic effects, and a theory that corrects for them would mean an 
improvement to the dynamical models for the motions of these 
satellites. 
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Table 6. Results for Titania 

Umbriel JPL 


Year 

Aq'* 

A5 

O-a* 

os 

ea* 

e6 

# 


(mas) 

(mas) 

(mas) 

positions 

1992 

-1 

-10 

27 

17 

54 

48 

69 

1993 

-8 

-32 

18 

22 

49 

54 

143 

1994 

43 

-58 

52 

54 

47 

48 

77 

1995 

34 

-20 

35 

38 

49 

50 

228 

1996 

-31 

-15 

67 

43 

55 

52 

285 

1997 

-22 

-20 

67 

25 

56 

47 

74 

1998 

3 

1 

33 

30 

49 

35 

29 

1999 

-57 

-1 

80 

57 

48 

49 

31 

2001 

-25 

-47 

34 

34 

39 

40 

65 

2002 

-11 

-6 

20 

23 

52 

50 

144 

2003 

-33 

-41 

61 

53 

65 

60 

118 

2004 

-28 

-48 

51 

41 

59 

56 

239 

2005 

38 

-65 

19 

15 

45 

52 

27 

2006 

-5 

-39 

70 

64 

56 

58 

82 

2008 

-10 

-62 

50 

56 

66 

57 

38 

2009 

19 

-115 

27 

31 

60 

59 

32 

2010 

-133 

43 

27 

8 

71 

67 

87 

2011 

-47 

-51 

46 

45 

59 

56 

820 

Columns: 

same as 

Table |3] 






Table 7. Results for Oberon 







Year 

Aof* A(5 
(mas) 

^ 6 

(mas) 

ea* es 
(mas) 

# 

positions 

1992 

2 

-8 

31 

16 

54 

48 

77 

1993 

-15 

-33 

20 

24 

48 

54 

119 

1994 

35 

-36 

25 

36 

50 

49 

64 

1995 

24 

-15 

40 

34 

49 

49 

230 

1996 

-36 

-9 

56 

33 

53 

49 

327 

1997 

-31 

-21 

100 

17 

56 

48 

93 

1998 

27 

7 

40 

23 

49 

36 

36 

1999 

-55 

-18 

66 

42 

48 

47 

36 

2001 

-21 

-39 

41 

31 

41 

43 

65 

2002 

-25 

-9 

19 

29 

53 

51 

154 

2003 

-39 

-17 

56 

47 

63 

61 

163 

2004 

-32 

-34 

43 

41 

58 

56 

295 

2005 

-29 

-56 

13 

17 

45 

53 

33 

2006 

-32 

-32 

58 

54 

55 

58 

142 

2008 

-21 

-60 

46 

50 

67 

58 

39 

2009 

55 

-88 

33 

38 

60 

59 

39 

2010 

-no 

42 

31 

12 

71 

68 

101 

2011 

-62 

-41 

44 

41 

57 

55 

915 


Columns: same as Tabled 
Table 8. Overall offsets, all satellites 



Umbriel JPL 



Fig. 7. Distribution of the offsets for Umbriel as a function of time. 


Titania JPL 



Titania JPL 



Fig. 8. Distribution of the offsets for Titania as a function of time. 


Satellite 

Aa* A6 
(mas) 

(mas) 

# 

positions 

Miranda 

-22 

-8 

96 

60 

584 

Ariel 

-30 

-21 

65 

40 

1710 

Umbriel 

-28 

-27 

62 

48 

1987 

Titania 

-25 

-35 

59 

48 

2588 

Oberon 

-35 

-26 

56 

42 

2928 


Columns: see Tabled 


5. Comparisons with other ephemerides 

We also compared our positions with those obtained from the 
recent planetary ephemeris INPOP13(0and from the dynamical 
model for the motion of the five main satellites of Uranus NOE- 
7-201^ This model is a more recent adjustment of the motions 
of the satellites than that given by Lainey (2008) 


^ See http://www.imcce.fr/inpop/ for the ephemeris and the respective 
scientific notes. 

® ftp.imcce.fr/pub/ephem/satel/NOE/URANUS/SPlCE/ 
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Oberon JPL 



Oberon JPL 



Fig. 9. Distribution of the offsets for Oberon as a function of time. 



Angle Sun-Observer-Uranus (degrees) 


Fig. 10. Offsets in right ascension (squares) and declination (triangles) 
as a function of the angle Sun-observer-Oberon. Each symbol is the 
average of 200 offsets. The horizontal line was set at ordinate 55 mas. 
Only data from Oberon were used here. 


In this section, only observed positions of Oberon are used 
for the comparisons. 


5.1. INPOP13C 


Table 10 lists the yearly offsets in the sense observed minus 
ephemeris (INPOP13c+ural 11) positions. We note that the right 
ascension between INPOP 13c and our observations agrees better 
than in DE432. However, the opposite is observed for the offsets 
in declination. The very last line of Table [T^ shows informa¬ 
tion that is equivalent to that in Table|^for Oberon, according to 
INPOP 13cH-ural 11, and can be assumed as a proxy to the offsets 
of the other satellites. 


B 

<1 


200 

Miranda + 

150 - Ariel X 


100 - 
50 
0 
-50 
-100 
-150 - 
-200 


+ 

U FI 


Umbriel 

Titania 


□ 


A M 


Oberon 


X + 


o^c^a^a^o^o^c^a^a^oooooooooo-^-^^ 

o^o^a^o^o^c^c^a^a^ooooooooooooo 

-H—i-H-H-^(M(N(N(N(N(N<N(N(N(N(N<N(N 


Date (year) 


to 

<1 


200 

Miranda + 

150 - Ariel X 


100 

50 

0 

-50 

-100 

-150 

-200 


Umbriel X 

Titania □ 


Oberon 


□ u 

+ 


s + 




-^fNfo•^^r)'sor-ooa^o—'(N<r)'^ir)'sC'r-ooavO-^(N 
o^o^a^c^o^c^c^a^a^oooooooooo-^-^^ 
o^c^a^G^o^c^o^a^a^ooooooooooooo 


Date (year) 

Fig. 11. Offsets in right ascension (upper panel) and declination (lower 
panel), as given by Tables to to all five satellites. Note that the 
offsets in right ascension for Miranda, measured in 2004 and 2006, fall 
outside the limits of the y-axis. 

Table 9. Overall residuals with respect to Oberon 


Satellite 

[Aof*] [Ad] 
(mas) 

(mas) 

# 

positions 

Miranda 

-2 

7 

69 

54 

584 

Ariel 

4 

1 

31 

28 

1710 

Umbriel 

3 

-6 

33 

30 

1987 

Titania 

5 

-6 

31 

30 

2588 


Columns 2 and 3 give the difference in the sense satellite off¬ 
set minus Oberon offset. Columns 4 and 5 give the respective 
standard deviations of these offsets. 


5.2. NOE-7-2013 

Table [m gives the yearly offsets in the sense observed minus 
ephemeris (DE432H-NOE-7-2-13) positions. The results are sim¬ 
ilar to those given in Table This is expected since, as men¬ 
tioned earlier in the text, the main source of differences between 
observed and ephemeris positions is the planetary ephemerides, 
which are responsible for the description of the motion of the 
barycenter of the Uranus system around the barycenter of the 
solar system. 

When we compare the offsets presented in Tables [T^ and [7] 
with those presented in Tables [TT| and [T) we note that the plane- 
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Fig. 12. Distribution of the offsets for Miranda as a function of the true 
anomaly. Solid dots: orbit view close to face-on. Note that there were 
not enough positions of Miranda with an orbit view close to the edge-on 
to appear in these plots. 


Fig. 14. Distribution of the offsets for Umbriel as a function of the true 
anomaly. Solid dots: orbit view close to face-on. Solid triangles: orbit 
view close to edge-on. 
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Fig. 13. Distribution of the offsets for Ariel as a function of the true 
anomaly. Solid dots: orbit view close to face-on. Solid triangles: orbit 
view close to edge-on. 


tary ephemerides clearly agree less well between each other than 
the models for the satellite motions (at least for the bodies stud¬ 
ied here). 


Fig. 15. Distribution of the offsets for Titania as a function of the true 
anomaly. Solid dots: orbit view close to face-on. Solid triangles: orbit 
view close to edge-on. 


6. Previous data from the Pico dos Dias 

Observatory for the main sateiiites of Uranus 

Raw images of the Uranus system, made at the Pico dos Dias ob¬ 
servatory before 1992, are no longer available so that VOS was 
the natural choice to provide an extension of this work to previ¬ 
ous epochs. 
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Fig. 16. Distribution of the offsets for Oberon as a function of the true 
anomaly. Solid dots: orbit view close to face-on. Solid triangles: orbit 
view close to edge-on. 


Table 10. Results for Oberon with INPOP 13c-l-ural 11 


Table 11. Results for Oberon with DE432-I-NOE-7-2013 


Year 

Aa* 

A5 

o-Q.* 

os 


£6 

# 


(mas) 

(mas) 

(mas) 

positions 

1992 

-2 

-9 

31 

16 

54 

48 

77 

1993 

-13 

-34 

20 

24 

48 

54 

119 

1994 

35 

-37 

25 

36 

50 

49 

64 

1995 

26 

-16 

40 

34 

49 

49 

230 

1996 

-34 

-9 

54 

33 

53 

49 

327 

1997 

-27 

-22 

97 

17 

56 

48 

93 

1998 

27 

7 

42 

23 

49 

36 

36 

1999 

-51 

-17 

66 

42 

48 

47 

36 

2001 

-19 

-40 

40 

30 

41 

43 

65 

2002 

-19 

-9 

18 

28 

53 

51 

154 

2003 

-39 

-19 

55 

47 

63 

61 

163 

2004 

-30 

-34 

45 

41 

58 

56 

295 

2005 

-33 

-56 

13 

17 

45 

53 

33 

2006 

-33 

-33 

59 

53 

55 

58 

142 

2008 

-28 

-61 

46 

50 

67 

58 

39 

2009 

55 

-90 

34 

37 

60 

59 

39 

2010 

-103 

44 

31 

12 

71 

68 

101 

2011 

-63 

-42 

45 

40 

57 

55 

915 

All years 

-34 

-27 

57 

42 

- 

- 

2928 


Columns: same as Table|3] 


Table 12. Overall residuals with respect to Oberon from V03 


Year 

Aa* 

^6 

a-a* 

0-6 

ea* 

es 

# 


(mas) 

(mas) 

(mas) 

positions 

1992 

-6 

-45 

31 

16 

54 

48 

77 

1993 

-20 

-70 

20 

24 

48 

54 

119 

1994 

35 

-72 

25 

37 

50 

49 

64 

1995 

28 

-51 

40 

34 

49 

49 

230 

1996 

-28 

-43 

56 

33 

53 

49 

327 

1997 

-19 

-54 

100 

17 

56 

48 

93 

1998 

43 

-24 

40 

23 

49 

36 

36 

1999 

-35 

-47 

66 

42 

48 

47 

36 

2001 

5 

-63 

41 

31 

41 

43 

65 

2002 

4 

-30 

19 

29 

53 

51 

154 

2003 

-7 

-35 

56 

47 

63 

61 

163 

2004 

2 

-49 

43 

41 

58 

56 

295 

2005 

8 

-68 

13 

17 

45 

53 

33 

2006 

6 

-41 

57 

54 

55 

58 

142 

2008 

19 

-63 

46 

50 

67 

58 

39 

2009 

95 

-89 

33 

38 

60 

59 

39 

2010 

-69 

44 

31 

12 

71 

68 

101 

2011 

-20 

-37 

44 

41 

57 

55 

915 

All years 

-8 

-42 

52 

42 

- 

- 

2928 


Columns: same as Table|3] 


Our aim here is not to change the values of the positions 
given by VOS. Instead, we simply applied two of the filters de¬ 
scribed in Sect. 4 to somewhat homogenise their position^ and 
ours. Procedures equivalent to the other filters were taken into 
consideration by these authors when preparing their data. We 
also provide a brief analysis of the respective filtered positions. 


^ Available at http://www.imcce.fr/hosted_sites/saimirror/butpomaf.htm. 


Photographic 

Satellite 

[Act*] [Ad] 
(mas) 

O'q'* (T s 

(mas) 

# positions 
before after 

Miranda 

-32 

-15 

164 

118 

315 

305 

Ariel 

-7 

-12 

66 

55 

363 

340 

Umbriel 

-11 

-4 

73 

60 

363 

335 

Titania 

-1 

-8 

52 

44 

363 

337 

CCD 

Miranda 

1 

6 

71 

72 

1130 

897 

Ariel 

4 

12 

39 

42 

1360 

1104 

Umbriel 

3 

10 

37 

51 

1350 

1113 

Titania 

-10 

17 

52 

61 

1365 

1214 


Filters described in steps 1 and 4, in this order (see Sect. 4), were 
applied separately to photographic and CCD data. Columns 1 to 
5 have the same meaning as those in Table Columns 6 and 
7 show the number of measurements before and after applying 
the second filter (based on the relative positions of Oberon with 
respect to those of the other satellites). Ephemeris positions were 
obtained from DE432-i-ural 11. 


In this context, we applied the filter based on the distance 
from the centre of Uranus (Perkin-Elmer telescope) and then the 
filter based on the relative positions of Oberon with respect to 
those of the other satellites. 

Eor a first verification, however, the data reported in VOS 
were separated into two sets: photographic (1988 and older) and 
CCD (after 1988) observations. The results are presented in the 
same way as for Table |^in Table [12] The filter related to the dis¬ 
tance from the centre of Uranus only eliminated Miranda data 
(62 positions) because the satellites were observed when their 
orbits were mostly face-on. 
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Table 13. Comparison between the positions reported by V03 and those 
from this work 


Original offsets in V03 

Satellite 

Aa* 

A6 

(Ta* 

0-6 

# positions 


(mas) 

(mas) 


Miranda 

-119 

-28 

114 

131 

326 

Ariel 

-128 

-38 

102 

123 

469 

Umbriel 

-120 

-46 

100 

117 

474 

Titania 

-136 

-44 

107 

115 

513 

Oberon 

-127 

-58 

106 

133 

551 

V03 minus this work 

Miranda 

-22 

17 

102 

133 

326 

Ariel 

-34 

21 

87 

120 

469 

Umbriel 

-33 

12 

90 

115 

474 

Titania 

-45 

17 

97 

116 

513 

Oberon 

-32 

0 

89 

130 

551 


Offsets from V03 (top table) and offsets from the difference in 
the sense positions from V03 minus (common) positions from 
this work (bottom table). 


< 
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I 





It is expected, in Table 12 that the standard deviations for 
the photographic measurements be larger than those made with 
CCD. In addition, since this work prohted from a better reference 
catalogue for astrometry, we consider that the efforts of VOS to 
determine relative angular measurements for the satellites were 
efficient, although the standard deviations in Table 12 are larger 
than those in Table 


In this context, we consider that both photographic and CCD 
positions are a homogeneous set of data, despite the much larger 
uncertainties of the Miranda photographic positions. From this 
point on, we therefore emphasise that the filter based on the rel¬ 
ative positions of Oberon with respect to those of the other satel¬ 
lites is applied without distinguishing between the detectors for 
the VOS data. This slightly changes the total number (CCD plus 
photographic) of positions for each satellite as given by the last 
column of Table ffel 

Table [TS] compares the original, hltered positions from VOS, 
with those obtained in this work. We note that VOS is closer to 
our results than to those of the ephemerides they used for com¬ 
parison (DE40S and GUST86). 

It is worth mentioning that VOS used the USNO-A2.0 
(Monet 1998) catalogue as reference for astrometry and applied 
corrections on it to minimise systematic effects on its positions. 
It is known from Assahn et al. (2001) that these effects are larger 
(in absolute value) than the offsets between VOS and this work 
(see Table [H- Therefore, the correction that VOS applied on 
the USNO-A2.0, as expected, brought the celestial frame origi¬ 
nally represented by that catalogue closer to that represented by 
the UCAC4. Most probably, smaller offsets and standard devia¬ 
tions would have been seen in Table pj] (VOS minus this work 
section) if proper motions were available to all reference stars 
used by VOS. These authors reported geocentric positions and 
we changed our positions accordingly for these comparisons. 

Figures[^to|^show the differences in the sense VOS minus 
this work as a function of time. The dimensions of the axes in 
these hgures facilitate the comparison with the similar axes in 
VOS. 


Fig. 17. Updated offsets of VOS for Miranda as a function of time. 
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Fig. 18. Updated offsets of VOS for Ariel as a function of time. 


7. Position cataiogues: Uranus and its main 
sateiiites 

Catalogues containing all the observed positions of the hve main 
satellites and Uranus, as well as the X and Y CCD coordinates 
of the observed satellites and reference stars, are only available 
at the Centre de Donnees astronomiques de Strasbourg (CDS). 
Table [T4| shows an extract of one such catalogue of positions for 
Oberon. Positions of Uranus (see TablepAlfor an extract) are also 
provided. Total counts are provided in Table [T6| 
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Fig. 19. Updated offsets of V03 for Umbriel as a function of time. 


Fig. 21. Updated offsets of VOS for Oberon as a function of time. 
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Fig. 20. Updated offsets of VOS for Titania as a function of time. 
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Since Uranus is saturated in most of the images, its posi¬ 
tions (see Table [T5] l were not measured. They were derived by 
applying the mean offset on the ephemeris position of Uranus as 
obtained from the satellites in the same image (except for Mi¬ 
randa). We stress that the positions of Uranus provided in this 
paper are not observed ones. 

We also note that ephemeris positions of Uranus and the off¬ 
sets for the satellites are both determined from DE432-i-uralll 
(this is also valid when we determine positions of Uranus from 
VOS data). Taking into consideration that the positions in Ta¬ 
ble were most frequently derived from the offsets of two to 


Table 14. Positions of Oberon - extract 


ID 

R.A. (ICRS) DEC. 

o-a* 


(UTC) 

(h m s) 

(O , 

(mas) 

2448782.68756100 

19 14 54.182 

-22 45 53.11 

57 

49 

2448782.70138981 

19 14 54.042 

-22 45 53.33 

56 

48 

2448782.73052454 

19 14 53.747 

-22 45 53.75 

65 

56 


Columns are Julian date (UTC) of the observation; observed 
ICRS right ascension; observed ICRS declination; uncertainty in 
right ascension - note the term cosd; uncertainty in declination. 
For the positions of Oberon and the other satellites from VOS, 
two extra columns provide the values to transfer their original 
geocentric positions to the topocenter. See text for details. 


four satellites, discrepancies in uralll (for instance, Figs.[T2|to 
[T6] | are attenuated and, therefore, the derived positions of Uranus 
can be regarded as ephemeris independent. 

Positions from VOS from before 1992 for Uranus and the 
main satellites are provided in separate catalogues (and sepa¬ 
rately from our data) with the same format as shown in Tables[T4| 
and[^ respectively. They are filtered by the procedure described 
in Sect. 6 and presented in their original form, that is, the right 
ascensions and declinations are exactly those provided by their 
paper. The application of the offsets given by the lower part of 
Table [TS] on VOS positions is a way to improve their results. We 
also note that two extra columns are provided in the tables, con¬ 
taining the data from VOS. They give the values Aacos^ and Ad, 
in units of milliarcseconds, that should be added to the right as¬ 
cension and declination, respectively, to transfer these positions 
from the geocenter to the topocenter. 

The observational data in VOS span from 1982 to 1998, so 
that all positions we provide after 1998 (5,827 for the hve satel¬ 
lites) are certainly unprecedented. 
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Table 15. Positions of Uranus - extract 


JD 

R.A. (ICRS) DEC. 

(Ta^ 

0-6 

sat. 

(UTC) 

(h m s) 

(O , 

(mas) 


2448782.68756100 

19 14 54.553 

-22 45 10.28 

57 

49 

o 

2448782.70138981 

19 14 54.431 

-22 45 10.51 

47 

42 

tauo 

2448782.73052454 

19 14 54.171 

-22 45 10.99 

58 

51 

to 


Columns are Julian date (UTC) of the observation; observed 
ICRS right ascension; observed ICRS declination; uncertainty 
in right ascension - note the term cosd; uncertainty in declina¬ 
tion; satellites that contributed to the determination of Uranus’ 
position; (a) - Ariel, (u) - Umbriel, (t) - Titania, (o) - Oberon. 
For the positions of Uranus from VOS, two extra columns pro¬ 
vide the values to transfer their original geocentric positions to 
the topocenter. See text for details. 

Table 16. Total counts for the catalogues 


Object # positions 

This work (from 1992) VOS (before 1992) 


Miranda 

584 

476 

Ariel 

1710 

474 

Umbriel 

1987 

462 

Titania 

2588 

558 

Oberon 

2928 

732 

Uranus 

3516 

732 


Number of positions for each satellite and Uranus. 


We stress that our positions are topocentric, whereas the orig¬ 
inal positions of VOS are geocentric. 


7.1. Uncertainties in position 

Uncertainties in position (see Tables [T4| and [T5]l are given by 

.6 — 


cr, 


icrl -I- al. 


( 1 ) 


For all satellites whose positions were determined in this 
work, cTa is the standard deviation in right ascension or declina¬ 
tion of the reference stars in the image containing the respective 
satellite and (Tt, is read from Cols. 4 (right ascension) or 5 (decli¬ 
nation) in Tablej^for the corresponding satellite. For Oberon, the 
values of S2 mas (right ascension) or 29 mas (declination) were 
used for cr^. They represent the mean values of the uncertainties 
in Table l^for Ariel, Umbriel, and Titania. For satellite positions 
from VOS, uncertainties were determined by taking cra=100 mas 
(right ascension and declination), and cr^ is read from Cols. 4 
(right ascension) or 5 (declination) in Table 12 for the corre¬ 
sponding satellite. For Oberon, the values of 43 mas (right ascen¬ 
sion) and 51 mas (declination) were used for crh. They represent 
the mean values of the uncertainties in Table [T^ (CCD section, 
where most of the positions are found) for Ariel, Umbriel, and 
Titania. The value of 100 mas for era is a tentative agreement be¬ 
tween Cols. 4 and 5 in Table [TS] (VOS minus this work section), 
in an attempt to reproduce the standard deviation of the reference 
stars in VOS. 

For all positions of Uranus determined in this work, cr^ is 
the standard deviation of the reference stars in right ascension or 
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declination in the respective image containing Uranus, and cr/, is 
the standard deviation of the satellite offsets in either right as¬ 
cension or declination, used to correct the ephemeris position of 
Uranus. When only one satellite was used, 32 mas (right ascen¬ 
sion) or 29 mas (declination) were used for cr^, see the reason¬ 
ing given above for these values. For positions of Uranus from 
VOS, uncertainties were determined by taking again cr„=100 mas 
(right ascension and declination) and crh as the standard devia¬ 
tion of the offsets of the satellites, in either right ascension or 
declination, used to correct the ephemeris position of Uranus. 
When only one satellite was used, the values of 43 mas (right 
ascension) or 51 mas (declination) were taken, see the reasoning 
given above for these values. 


8. Conclusions and comments 


We determined accurate positions of the five main satellites of 
Uranus: Miranda, Ariel, Umbriel, Titania, and Oberon for the 
time span 1992 - 2011 from observations made at the Pico dos 
Dias Observatory. Positions of Uranus, derived from those of the 
satellites, were also determined. The standard deviation of the 
offsets is typically between 40 mas and 65 mas. 

Our positions contribute to the general knowledge of the 
physics of the Uranus system (see Jacobson 2014) to improve 
modem planetary ephemerides, such as DE432 and INPOPlSc, 
as well as modern models of the satellite motions, such as ural 11 
andNOE-7-2013. 

A number of observations made far from the opposition are 
certainly useful to improve the determination of the heliocentric 
distance of Uranus. 

The positions from VOS can be used to provide an extension 
to previous epochs of our data. The application of the offsets pro¬ 
vided by the lower section of Table 13 on their original positions 
is a way to improve them. 
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